HD-THEP-92-50 

TPI-MINN-92/62-T 

TUM-T31-32/92 

Remarks on the Heavy-Quark Symmetry 
in the Complex Plane 

Patricia Ball 

Physik- Department, TU Miinchen, D-8O46 Garching, FRG 

H.G. Dosch 

Institut fiir Theoretische Physik, Universitdt Heidelberg, 
D-6900 Heidelberg, FRG 

M.A. Shifman 

Theoretical Physics Institute, University of Minnesota, 
Minneapolis, MN 55455, USA 

February 1, 2008 



Abstract 

Hadronic systems built from a heavy quark and a cloud of light quarks and 
massless gluons posess the Isgur-Wise symmetry resulting in certain relations 
between various form factors. These relations can not be valid in the entire 
complex plane. We identify contributions which unavoidably violate strongly 
the symmetry in the time- like domain near thresholds. A possible connection 
with the heavy-quark phenomenology is discussed in brief. 



submitted to Physical Review D 



1 Introduction 



Strong interactions in quark-antiquark systems become simpler if the mass of one of 
the quarks goes to infinity In particular, the transition amplitudes of the type 
{Qq) (Q'q), where Q is a generic notation for a heavy quark with mass tuq, which 
will be used throughout the paper, and q denotes a light (massless) quark, induced 
by the vector and axial vector currents Q'-^^Q and Q'-^^'^^Q, respectively, are related 
by a symmetry which takes place even if mq ^ rriQi, but both masses, mq and mg/, 
are parametrically larger than Aqcd- This symmetry is called the Isgur-Wise (IW) 
symmetry 0. Moreover, the absolute normalization of the form factors in the small 
velocity limit is fixed |^.[] 

For finite quark masses the symmetry is, obviously, broken. The parameter which 
governs the breaking of the IW symmetry can be represented by some positive power 
of the quantity 

/If 2 _ nf2 

M| ■ 

Here B* and B denote generic vector and pseudoscalar mesons (ground states) with 
the quark content (Qq). The quadratic mass difference 

5^ = Ml, - Ml (2) 

stays constant in the limit niq oo; therefore, one might think that in this limit 
the heavy quark symmetry (essentially, independence of the strong interactions of 
the spin orientation of the heavy quark) becomes perfect. While in a large range of 
momentum transfer this should be the case, there exist kinematic domains where e 
does not actually measure the strength of the symmetry breaking. This happens, 
in particular, when the form factors are considered in the time-like region, near 
the corresponding thresholds. Under these circumstances the genuine symmetry 
breaking parameter is 5'^ /{q^ — 4M^) rather than e , and relations stemming from 
the IW symmetry do not hold. As a matter of fact, in certain instances the breaking 
is even parametrically large. This simple observation will be quantitatively worked 
out below. 

This phenomenon is not unique, of course. Similar behavior is well-known, say, 
for the isotopic invariance of strong interactions. Amplitudes which are, generally 
speaking, believed to be isotopically invariant can exhibit very strong deviations from 
the symmetry predictions provided that a typical energy scale in the process at hand 
is not large compared to the mass difference m^ — m„. A classical example is the DD 
form factor in the vicinity of 7/^(3.77). The imaginary part of this form factor in the 
resonance is strongly different in the channels D^D^ and D^D^ due to the fact that 
the energy relase is about 40 MeV - quite comparable to 2M{D^) — 2M{D^) ^ 10 



^Historically the observation ||] that the transition form factors are predictable in the small 
velocity limit for uiq ^ itiqi served as an initial impetus for the introduction of the IW symmetry 
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MeV. It is quite typical that the symmetry breaking effects, although large, are 
calculable in this case. 

In this note we analyse, in the same spirit, the peculiarities of the IW symme- 
try breaking in the time-like region, mostly due to anomalous thresholds @, A 
remark on possible phenomenological implications in the so-called molecular quarko- 
nium (analog of the molecular charmonium P, 0) is presented. For definiteness we 
concentrate on the diagonal vector current 

= Qlf^Q, (3) 

although the assertions made below are of a general nature and are applicable, in 
principle, to other currents - axial, non-diagonal, etc. For the vector current one 
can consider three transition amplitudes: 

< B\J^\B >, < B\J^\B* > and < B*\J^\B* > . (4) 

Taking into account current conservation one concludes that the first and the second 
transitions are described by one form factor each, 

< B{p')\J,\Bip) >=F+P„ 

< B{p')\J^\B*{p) >= -ifSf^^fs^p'^Pfse^, (5) 
while the third amplitude is represented by four form factors, 

< B*ip')\J,\B*{p) >= F^{ee')P, + F,[e,{e'P) + e'^{eP)] 
+ Fs^-^^^P, + F,[e^ie'P) - e'^Pf^^^^j^^, (6) 

where 

P = p' + p, q = p' -P, (7) 

is the B* mass, e and e' are the polarization vectors of the initial and final i?*'s. 
All six form factors, F+, /, Fi, F4, are functions of g^. 

The heavy quark symmetry reduces the number of independent functions in the 
general decomposition (5) and (6) to one, in the limit rriQ —>■ 00 one obtains : 

j^<B\J^\B>=av' + v)„ 

1 

— < B\J^\B* >= -i^e^apWa^f^^^i^ 

^ < BV,\B* >= a-{^e'){v' + v), + [e,(e'K + v)) + e'^{e{v' + v))]]. (8) 

Here = p^/M, and the function ^ depends in the IW limit only on the product y 
of four-velocities, 

y = v'v. (9) 
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Moreover, ^(y = 1) = 1 0. Let us draw the reader's attention to the fact that 
and F4 are predicted to vanish. 

Any deviation from these relations will signal a violation of the IW symmetry 0. 
An obvious violation occurs above the threshold of the BB production but below 
BB*. Indeed, in this domain all form factors have imaginary parts associated with 
the normal thresholds due to BB. On the other hand, there is no contribution to 
the imaginary part from BB* and B*B*. It is quite clear that in the pseudoscalar 
B meson the spin of the heavy quark Q is rigidly correlated with that of the light 
cloud. Hence, the spin independence of the heavy quark interaction - raison d'etre 
of the IW symmetry and the origin of eq. (8) - is totally lost. In particular, the 
"forbidden" function F3 appears. Below we will demonstrate this assertion in a 
less trivial context of anomalous thresholds. The existence of these thresholds is 
due to pion exchange. A specific feature which makes their analysis interesting 
is the fact that they can start parametrically much below the normal thresholds, 
depending on the interplay between 6 and m^, the pion mass. Certainly, in the 
real world both quantitites are such as they are. It is instructive, however, to play 
with m-Tr/^ treating this ratio as a free parameter. One can easily change by 
adjusting the light quark mass terms in the QCD lagrangian. We will investigate 
different regimes allowed for the anomalous thresholds in Sect. 2. Among other 
things it will be seen that under certain conditions the anomalous threshold starts 
at a value of which is independent of the heavy quark mass. The corresponding 
contribution, obviously, has no smooth IW limit in the sense that its dependence 
on the momentum transfer does not reduce to a y dependence, as prescribed by 
eq. (9). This particular contribution is not leading at tuq 00. There are other 
regimes, however, (which also defy the heavy quark symmetry) where we find an 
enhancement by positive powers of rriQ in the near-threshold domain. 

The issue of molecular quarkonium (bound states of the type B*B*) and its 
possible relation to the violations of the IW symmetry is discussed in Sect. 3. 



2 Pion Exchange and Anomalous Thresholds 

First of all it is instructive to notice that the form factor of B is free from anomalous 
thresholds while that of the B* meson is not. Indeed, let us consider the triangle 
graphs depicted in Figs, la, b. Both graphs can lead to singularities on the physical 
sheet below the normal thresholds the positions of which depend on the external 
masses. The existence of such anomalous singularities has been first realized by 
Karplus, Sommerfield and Wichman Q. From the standard Landau theory |jlO| it 



is not difficult to find the conditions under which the anomalous thresholds may 
occur. 



^The above reduction of six form factors to one for non- relativistic heavy quarks, i.e. at \v\ <C 1, 
is known in the literature for more than fifteen years M, see also Chapter 4 in the review paper 
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Inspecting the dual diagrams associated with Figs, la, b we conclude the fol- 
lowing: 

(i) The graph la has an anomalous threshold provided that 



The beginning of the anomalous singularity is at 



(10) 



BB 



4M^ 



:iii 



To avoid numerous sub- and superscripts we have introduced the notations t = q"^, 
fi = ttIt^, M = Mb and = Mb*- The position of the anomalous threshold varies 
from zero at the lower boundary of the interval (10) to the normal threshold 4M^ 
at the upper boundary. 

(ii) The condition on /i guaranteeing the existence of an anomalous threshold in 
graph lb is: 

< ^ < . (12) 



M + M, 



1 + (M/M,) 



The lower and upper boundaries in eq. (12) correspond to the position of the anoma- 
lous threshold at M^M, + M) and {M + M^f, respectively. 

The general expression for the position of the anomalous threshold in this case 

is 



M^ + Ml + -{5^-li'') + 



AMI 



-l){^^''M^ 



,5^ -fx' 



nl/2 



(13) 



A remarkable feature of eq. (11) is the fact the the position of the anomalous 
threshold, t^^ , needs not to scale as in the limit mg oo; if /z is chosen 
to lie close to — M other regimes are quite possible. For instance, if /i = 
(M* — M)(l + const. M'"^) we find that t^^ scales like M°. The momentum transfer 
dependence then will be expressed by a function of (M^/t^^)(l — y) rather than a 
function of y, and the IW scaling will be obviously violated. Of course, the above 
regime is exotic since it requires fine-tuning of the pion mass. One can always fix 
the pion mass and then proceed to infinitely heavy quarks. This will push to to 4M^ 
and will restore the y dependence. 

The position of the anomalous threshold in the complex plane, although interest- 
ing by itself, says nothing about the relative weight of the anomalous singularities. 
A priori it is conceivable that they totally decouple at rriQ oo. Whether or not 
they actually decouple depends on the behavior of the pion constant. 

A brief reflection shows that the leading M dependence of the pion coupling to 
the heavy mesons should be such that after proceeding to non-relativistic normal- 
ization M should disapper altogether [0, |ll|. A relation between the B*B*Tr and 
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B*Bt[ vertices can be readily obtained by combining the heavy quark and chiral 



symmetries [12, 13, O]. In the relativistic normahzation one gets 



B^B-k: ^M{ek) +0{M'>), (14) 



B*B*7r: -i—e^p^s{p + p')^kpe'^e5, (15) 



r: 

where p , p' and k = p — p' are the momenta of the initial and final heavy measons 
and the pion, is a dimensionless coupling constant of order 1 ||12| , [T3| , [T^ , and 
is the pion decay constant. It is easy to check, for instance, that with eq. (p!^ ) the 
D* — > D-K width does not contain the heavy meson mass, in full accordance with 
the arguments of refs. ITT . 



We pause here to make an important remark. The fact that the B*B7i vertex is 
proportional to M (in the relativistic normalization) is a consequence of the standard 



assumptions of the effective heavy quark theory KB, 16|. Under certain kinematic 



conditions when the IW symmetry is broken in the sense discussed in this note. 



the BB*7i vertex may deviate from the standard scaling law exhibited in eq. (0). 
In other words, this vertex may be proportial to another power of M for trivial 
kinematic reasons which will become clear shortly. For the time being let us stick. 



however, to eq. (0), with g = 0(1). 

It is important that the discontinuity of the amplitude at the anomalous cut 
is determined by the triangular graphs of Fig. 1 with all three internal particles 
on mass shell 0, |T^. Therefore, the pion coupling constant which enters all our 
formula is the on-mass-shell coupling. 

Since our purpose is mostly illustrative it is reasonable to limit our analysis to the 
anomalous singularities of graph la .The singularity of this graph can, in principle, 
be arbirary close to t = while that of Fig. lb necessarily lies higher than 2M^ (see 
eq. (13)) and, moreover, can be pushed to the unphysical sheet provided that /i is 
chosen larger than S{1 + (M/M^))-^/^. 

The contribution of the triangle graph depicted on Fig. la can be written as 
follows: 

, _1..,2 / 1 {ke)ike'){p + p' + 2k), 

^ J (27r)4A;2-/i2((p+A;)2-M2)((p' + A;)2-M2)' ^ ^ 

where according to eq. (14): 

S-f. (17) 

In 

The anomalous cuts can be easily evaluated using Cutkosky [|1^ rules, yielding 
the imaginary parts pi(t) of the form factors Fi(t), 

= 3 Vt(4M2 _ ,)5/2 ' (18) 
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^ 2 ml^^' -6'- ^^\25^ + + t)){25' + 
^'^^ ^ Vt(4M2-t)5/2 ' ^ > 

p,{t) = -2M,2^pi(t), (20) 

Pa = 0. (21) 

It is quite obvious that the above expressions badly violate the heavy-quark 
symmetry, since pi ^ —p2 and, moreover, ps ^ 0. The above equations should be 
supplemented by the statement that the anomalous part does not vanish in the limit 
M — i> oo. It is convenient to carry out the analysis of the absolute normalization 
separately in two distinct cases: 

(i) The case in which 6 and p are fixed and thus the anomalous threshold is 
near the normal one. This is referred to as the relativistic case; 

(ii) the case in which 6 stays fixed but p decreases at least as and thus 
the anomalous threshold is far below the normal one. This is referred to as the 
non-relativistic case. 

It should be noted that in the case considered here the anomalous singularities are 
due to hadronic intermediate states and not to quarks, as in the case of formfactors 
of heavy quarkonia, considered in refs. []TB|, . 



2.1 Relativistic case 

If 6 and p are fixed the anomalous threshold stays, according to eq. (11), in the 
vicinity of the normal threshold, i.e. the relevant variable is m = 4M^ — t. Keeping 
only the leading in M terms we get 

. . _ ,2 (2^^ + 2p^-n)(^^ + 2^V + p^(p^-n)) 
P^W - 9 (23) 



P3[U) = -9 



~2 Mt 



2MI'^ u7/2 

{pV - 3m((5^ + 4 (5V + 3p') + 10(5" + i?f} (24) 

The anomalous cut in the u plane extends from (5^+p^)^/ p^ to 45^. If m = 0(5^) 
and g = 0{M) the discontinuities are large rather than small, 

Pi(m), P2{u) (X M, P3(m) (X M\ (25) 

leading to a huge violation of the heavy-quark symmetry in the near-threshold do- 
main (by the near-threshold domain we mean an interval of t centered at 4M^ whose 
length is 0(5^).) Furthermore, it is quite trivial to do the dispersion integral with 
the discontinuities given in eqs. (22-24). With 

^ 1 ^ p„ 

71 J to t — t 
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one finds outside the near-threshod domain (at t = 0(M^)) 



^an _ pan _ pan _ (^27) 



Specifically we obtain 



F^-(f -0)- A^l^J^ - 2 (^^-/^TM. 



The expression (^8]), taken at its face value, would signal a parametrically large 
symmetry violation even at t = 0. We hasten to add that the anomalous terms we 
have calculated should be (better to say, are expected to be) cancelled by normal 
contributions of the hadronic graphs if one considers the form factors Fi outside the 
near-threshold domain. 

It is instructive to check how the above cancellation works by directly computing 
the graphs of Fig. 1 (this is not a realistic computation of the form factors for many 
reasons, of course, but just an exercise allowing one to see the restoration of the 
symmetry). If t is not especially close to the threshold, the diagram la is not singled 
out. One must consider all possible triangle graphs, both for B* and B, with the 
intermediate states containing BB, BB* and B*B* . Adding them together, we 
observe that all symmetry-violating terms with positive powers of M cancel, if we 
start with formfactors satisfying eq. (8). 

2.2 Non-relativistic limit 

Let us consider another limiting case, when /i is only slightly larger than — M, 
so that fi scales as M~^. Then the anomalous threshold begins far below the normal 
one, i.e. to < AM^. 
We start with 

= -^:r- (29) 
^ 2M7 ^ ^ 

where 7 is a fixed number close to 1 but slightly less than one. Then 

to ^ (1 - 7^)4Ml (30) 

Neglecting t compared to we get for the anomalous discontinuities near the point 
t = to the leading terms: 
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The anomalous contribution depends on t through t/[4(l — 7^)M^], and the form 
factors scale in the following way 

pan _ pan _ ^-2^ pan _ ^-4^ (^34^, 

In other words, the anomalous contribution decouples in the limit mq oo. 

Of more interest is the case when 7 — 1 for M 00. Specifically, let us consider 

/i=(M,-M)(l + ^), (35) 

where a; is a parameter of order 5^. Then ~ 8a; = 0(M°). Under this choice - 
when M + /i — M^, = E = 0{M~^) - it is reasonable to turn to a model according 
to which B* is a non-relativistic bound state of Bit, or, at least, such a four-quark 
(molecular) component is present in B* with a certain prabability. 

To make the situation more graphic, assume, at first, that B* completely reduces 
to a loosely bound system of Btt, analogous to the deuteron. Clearly the spin sym- 
metry between B* and B is maximally violated in this case, since the pseudoscalar 
meson, the would-be partner of the B*, does not look like this bound B* at all. 
Moreover, in this case one would expect that the heavy-quark symmetry is violated 
not only in the near-threshold domain, but everywhere in the complex plane. The 
question is how this comes out formally. 

The answer to this question is rather obvious. If B* is like the deuteron, its 
coupling constant g is rigidly fixed in terms of the binding energy E (see e.g. 



for an S-wave bound state) and turns out to be much larger than that given in eq. 
(17). Before we proceed to derive the relation between E and g let us quote the 
expressions for the anomalous cuts in this limit. 

If we had substituted g = 0{M) we would have got at t = 

Fr(0) ~ l/M^ F2""(0) ~ F^'^iO) ~ 1/M, (39) 

Now, let us derive the actual scaling law for g in this scenario of a loosely bound 
state. Consider the non-relativistic contribution of diagramm la , by doing the 
corresponding manipulations in eq. (16). We put the pion on mass shell, i.e. take 
only the imaginary part of its propagator, which just singles out the anomalous 
contribution. We then replace the zero components of the occuring Lorentz-vectors 
by their corresponding nonrelativistic expressions: 

Po = M,+ p V(2M,) , ko = fi + P/{2fi) (40) 



8 



The pion is nonrelativistic, since the binding energy and hence the virtual momenta 
are small as compared to the pion mass according to (34). It is most convenient to 
work in the Breit-frame, where p = —ff = q/2 . We also neglect the pion mass /i 
against the heavy mass M^, . This yields for the zero component of the current: 

nonrel _ ^9^ f ^3j. ij ' ^ ' /i/(2M,)g ■ e){^^ k - fi/{2M,)^ ■ q) ^^^^ 



16M,(27r)3 J {k- /i/(2M,)g)2 + a'){k + ^/(2M,)g)2 + 

where a = y/2fiE with the binding energy 

E = M + fi-M^ (42) 

We compare this expression to the form factor of the S*-meson, which is as- 
sumed to be a bound state of a 5-meson and a pion. We evaluate it in the impuls 
approximation with the momentum space wave function 

^{k) = NJ{k/a)k/{k'^ + a^) (43) 

with k = \k\ . The unphysical pole at k"^ = —o? reflects the exponential tail of the 
pion cloud. 

For the normalization constant we obtain 

W' = 4^"/(P^^?/' (44) 

i.e. N\ ^ C/a. The impuls approximation for the zero component of the vector 
formfactor yields: 

A*-"- = ^ y (i'k^ ■ t/j{k - /i/(2M,)g) rik + i^/{2M,)q) ■ e (45) 

A comparison of ( ^T]) and (|45| ) shows, that the function f{k/a) in ( ^3] ) corresponds 
to some momentum cutoff in (|4l|) , and that the two expressions coincide, if the 
coupling constant is given by 

g = (46) 

where the constant C depends on the momentum cutoff. Notice that the relation 
between g and E is different from that discussed in |2^ due to the fact that the 
resonance we consider is P-wave. 

Thus, for fi scaling like (^5]) we get 

f oc Ml (47) 

It is rather straightforward to check that the form factors within this scenario will 
have nothing to do with eq. (8), as expected. 
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A more realistic scenario will be to say that there is an admixture of the Bn 
bound state in B*. If the amplitude of this component scales as M~^/^, we have 
g = 0(M) and return back to eq. (PU]). 

The scale for the t-dependence of the anomalous contributions is now however 
t^^, which does not scale with and hence we have for that case a -not very 
spectacular - violation of IW symmetry, which might be remarked as a spike near 
t = in the derivatives of Fi{t). 

The question "what is the actual admixture of the molecular Bit state in B*7" 
is a dynamical one, and cannot be solved on the basis of essentially kinematical 
arguments presented above. The answer depends on the intricacies of the large 
distance dynamics and can be a non-trivial function of E. While intuitively it is 
clear that the overlap is less than 1, it needs not be as small as M~^/^. Then g will 
not scale as prescribed by eq. (17). If the overlap is larger than 0(M~^/^), the spin 
symmetry violation in the entire complex plane will survive in the limit M —>■ oo. 

We should mention that in the chiral limit n = there is no anomalous threshold 
on the physical sheet. In that case the B* is unstable, but its width vanishes as 

3 Molecular Quarkonium 

We would like to remind that the well-forgotten issue of the molecular quarkonium 
- loosely bound states in the systems BB^ BB* + B*B^ B*B* is relevant to the 
discussion of the heavy-quark symmetry in the complex plane. If the long-range 
forces binding these "molecules" are independent of the heavy quark spins 0], i.e. 
eq. (8) holds, it is easy to check that the ratio of yields in these channels is 1:4:7 
0. This is the ratio stemming from eq. (8) after squaring these expressions and 
averaging over the spacial orientations of the momenta. The question is how the 
symmetry is realized, if at all, in the "molecular" domain. (The "molecules" can be 
either bound states or resonances above threshold.) 

The answer depends on the relation between the spin splitting 6"^/ (2M) and the 
molecular- level splittings. For the lowest molecular levels the latter are expected to 
be parametrically larger than 5^/(2M). Then the spin symmetry will be reflected 
in an (approximate) triple degeneracy of the levels. On average (summing over the 
triplet of the degenerate levels) all symmetry relations, including 1:4:7 for the ratio 
of yields, will be fuUfilled, but at exactly the position of each individual resonance 
they are maximally violated. On the other hand, if there are molecular states 
very close to the threshold, Mres - 2M, < 5^/{2M), it can well happen that such 
resonances exist, say, only in the B*B* channel and are absent in BB (as is the case 
in charmonium). Then the form factors near the resonances will not be symmetric. 
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4 Concluding Remarks 



We have investigated the imphcations of anomalous thresholds of heavy meson 
formfactors due to pion exchange. For the actual values of the pion mass and 
the quadratic mass splitting of the pseudoscalar and vector heavy mesons we have 
found violations of the IW symmetry near the production threshold of pairs of these 
mesons. These violations are not only the trivial ones, due to the finite quadratic 
mass splitting, but singularities of formfactors, which are predicted to vanish become 
parametrically large, i.e. increase with the heavy quark mass mq. The relation of 
these violations of the IW symmetry near t — Atuq with the molecular quarkonium 
is stressed. 

If the pion mass is treated as free parameter, anomalous thresholds might ap- 
proach the region t — Q arbitrarily close. The quantitative importance of such 
"low-lying" singularities depends on the degree of admixture of the molecular Bn 
state in the i?*-meson. 
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Figures 




Figure 1. Diagrams witli pion excliange contributing to the B* formfactor. 
Note that there is also a diagram corresponding to Fig. lb with the internal B- and 
S*-lines interchanged. 
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